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Abstract. The plasma dispersion function is computed for a homogeneous 
isotropic plasma in which the particle velocities are distributed according 
to a Kappa distribution. An ordinary differential equation is derived for 
the plasma dispersion function and it is shown that the solution can be 
written in terms of Gauss’ hypergeo metric function in the form Z K ($) = 
iK~ 3 / 2 (K—l/2)F(l, 2 k, K+l,t), where t = (l+i£/-\//c)/2. Using the extensive 
theory of the hypergeometric function, various mathematical properties of the 
plasma dispersion function are derived including symmetry relations, series 
expansions, integral representations, and closed form expressions for integer 
and half-integer values of k. 


1 Introduction 

Particle distribution functions with power law tails are frequently ob- 
served in space plasmas throughout the solar system. In practice, they are 
often modeled using Kappa distributions (defined below). The existence of 
such nonthermal, non-Maxwellian distribution functions can have profound 
effects on wave propagation and transport processes in these plasmas; effects 
that are of significant interest in space physics. The kinetic theory of such 
processes depends in an essential way on the plasma dielectric function and, 
therefore, on the plasma dispersion function. 

The purpose of this paper is to derive the plasma dispersion function 
for the Kappa distribution in a new way and to investigate its mathemat- 
ical properties. This function was derived independently by Summers and 
Thorne 1 and by Chateau and Meyer- Vernet 2 for the special case where k = 
1, 2, 3, . . . It was later derived in the more general case k > 0 by Mace and 
Hellberg. 3 These authors all used the technique of contour integration. The 
approach taken here is different, being based instead on differential equations. 


In order to be as consistent as possible with the formulation of the classical 
plasma dispersion function, the definition given here is slightly different from 
that of Summers and Thorne 1 and Mace and Hellberg. 3 To clarify the rea- 
sons for doing this, the definition of the plasma dispersion function is briefly 
derived from first principles. The resulting dispersion relation for plasma 
waves, equation (29), is somewhat simplified in the present formulation. 
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Section 2 gives a concise list of the main results of this paper. The 
definition of the Kappa distribution is given in Section 3, and the definition 
of the plasma dispersion function is given in Section 4. In Section 5, the 
differential equation for the plasma dispersion function is derived and the 
solution is expressed in terms of Gauss’ hypergeometric function. Symmetry 
relations and series expansions are derived in Sections 6 and 7. Integral 
representations are derived in Sections 8 and 9. Closed form expressions 
for integer and half-integer values of k are derived in Sections 9 and 10, 
respectively. 


2 Summary of Results 

The important properties of the plasma dispersion function are listed here 
for easy reference. These properties are investigated in detail in the sections 
that follow. Throughout this paper the power function z a , a G M, is defined 
by its principal branch, that is, z a = r a exp(ia6), where z = rexp(i$), r > 0, 
and —7 r < 6 < tt. 

Plasma Dispersion Function 


Zn{C) 


D p 

k Z/2 


1 , 2/y k T 




K > 0 . ( 1 ) 


If k — n is an integer, this function has a single pole of order n at the point 
£ = —iy/P. Otherwise, it is analytic with a branch cut from £ = — 1 \[k to 
—ioc along the negative imaginary axis. 


Behavior Near £ = —i\fn 


Z K (0 ~ 


m 1/2 T(/c) / 

2- 1 v ^r(K- !) V 




( 2 ) 


Symmetry Relations 


z K {-C) = -[z K {0X- ( 3 ) 

2 r(/c) / £ 2 V* 

Kp 2 T(/c - 1/2) V 


Z K {0 + Z k {-£) 


K 


( 4 ) 



J.J. Podesta — Plasma Dispersion Function for the Kappa Distribution 3 


Integral Representations 
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Differential Equation 
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Differentiation Formula 
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Closed Form Expression for n — n — 1, 2, 3 . . . 
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where z — £/ \Jn. 

Closed Form Expression for n = n + | 
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where z = £/ y/H and n— 1, 2, 3, . . . 

Series Expansion for \z\ < 1 



where z — £/ y r if. 

Series Expansion for \z\ > 1 


^«(0 = yS T{h-±) { 1 + i tan (^) s g n [ Re (")] } (* 2 ) K ( 1 + 



where z = £/\/k, k ^ and the function sgn[Re(z)] is defined by 

equation (81). 


3 Kappa Distribution 


In the astrophysics and space physics literature, the isotropic Kappa dis- 
tribution in three dimensions is defined by 


v 


2 \-( K +!) 


/k(v) = A k H 2 » K> !/ 2 > 




(14) 


where v = ( v x ,v y ,v z ) is the velocity vector, v = |v|, and Vq is some charac- 
teristic velocity. The normalization constant is 


A K — 


T(k + 1) 


r(«) 


(tTKVq) 3 / 2 T(k— l /2) (tTVq) 3 / 2 ^ r(/C— y 2 ) 

and is chosen such that 


(15) 


4vr / f K (v)v 2 dv = 1. 


(16) 
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The condition k > 1/2 is necessary for this integral to converge. The Kappa 
distribution is closely related to the Beta distribution. In fact, the Beta func- 
tion arises when computing the normalization factor A K and the statistical 
moments. In the limit as k — > oo, the Kappa distribution approaches the 
well known Maxwell distribution e~ v /v °. 


The velocity moments of the Kappa distribution are given by 

r°° 9 ri w+3 iriv _ »+i \ 

(V}=4W Mv)v" +2 dv = 1 2 1 , 2 f (17) 

Jo r(« — 2) 


where n is an integer and 0 < n < 2(k — 1). For an arbitrary real power 
n > 0, this integral is finite if and only if n < 2k — 1. In applications, the 
equivalent temperature of an isotropic Kappa distribution is defined by the 
relation 


\m <„ 2 > = \k B T 


(18) 


which implies 


knT = -mvn 


K 


K 


3/2 


(19) 


where ks is Boltzmann’s constant. The second moment (v 2 ) is finite if and 
only if k > 3/2. 


4 Plasma Dispersion Function 


In the Vlasov theory of plasma waves 4 the propagation of electrostatic 
waves (having E j| k) in a homogeneous isotropic plasma is governed by the 
dispersion relation ^(k, s) = 0, where 


£l{ k, s) 



dF o0 /du 

u-is/\k\ 


( 20 ) 


is the longitudinal dielectric function and SI units are used throughout. The 
index a denotes the different particle species, c o pa is the plasma frequency 
for species a, k is the Fourier wavevector with magnitude k = |k|, s is the 
Laplace transform variable, and F a0 is related to the particle distribution 
function f a0 by 


k v 

~w 


drv. 


Faoiu ) 


fao(v)S(u 


( 21 ) 
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The distribution function / a0 (v) is normalized to unity. 

Spatial isotropy implies /o(v) = /o(u 2 ), where v = |v|. After a rotation of 
coordinates which aligns the v z direction with k, the last equation becomes 


FaoM = 



fao(vl + V l + u 2 ) dv x dv y 


= 2vr / fao{v 2 )vdv 
J\u\ 


and, therefore, 


dF ( 


0:0 


du 


= -2nuf a0 (u 2 ). 


(22) 

(23) 

(24) 


It is usually the case that f a 0 has the functional form f a o( v / v ao ), where v a0 
is some characteristic velocity; consequently, F a0 = F a0 (u/v a0 ). Changing to 
normalized velocity variables: 


u dF 1 dF 

u = — , F(u') = v a0 F(u'), — = — — , 

v a0 du r)r "’ 


v'io du' 


the dispersion relation (20) can be written 




\k\v M , 


(25) 


(26) 


where, after dropping the primes on u, 
3 dF a0 /du 


Z\a,0 = 


and, by definition, 




du = 


F a0 (u) d 

du — — 


F a o(u ) 


Z{a,i) = 


1-00 (« - 0 2 
00 F a0 ( u ) 


df J u-f 


'-oo u ~£ 


du, Im(£) > 0, 


du (27) 


(28) 


is the plasma dispersion function. It is defined for Im(£) < 0 by analytic 
continuation. Introducing the complex frequency u = is which is consistent 
with Fourier modes of the form e h k - x -wt), ^j ie pi asma dispersion relation takes 
the standard form 5 


i-E 


U pa ^ ^ 


k 2 u 2 a0 


Dk\v a0 ) 


= 0. 


(29) 
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If f a o is a Maxwell distribution 

fao(v) = 


-D/vl 


l TTVi 


;) 3/2 


(30) 


then the evaluation of (23) and (28) yields the well known plasma dispersion 
function 6 

1 f°° e - “ 2 

I jdu, lm(0 > 0. (31) 


l-oo u ~£ 

If fao is a Kappa distribution (14), then the evaluation of (23) and (28) yields 

1 r(/c) 1 


Fao(u) = 


and 


5.(0 = 


y/nifvo r(/c— |) (1 +u 2 /kVq) k 

T(k) f°° 1 


KFk T(k — \) oo (u -0(1 + u 2 /k)' 


du. 


(32) 


(33) 


This is the dispersion function for the isotropic Kappa distribution. It gen- 
eralizes the well known result (31) for the isotropic Maxwell distribution. 
It is different from the definition of Summers and Thorne 1 which was also 
adopted by Mace and Hellberg. 3 These authors employ an exponent k + 1 
rather than k in the integrand. Consequently, the function Z K { 0 of Summers 
and Thorne is more closely related to the derivative Z' K (f) of the function 
defined here. 


5 Representation by Gauss’ Hypergeometric 
Function 


For purposes of mathematical manipulation, it is convenient to let x = 
u/y/if in equation (33) and write 


ZM) 


-i_ r M r — —dx 

i/sr(K-!)F(i + i ! )' (x-z) 


Im(^) > 0, (34) 


where z = £/>/&■ Therefore, the calculation of the plasma dispersion function 
reduces to the evaluation of the integral 


I 



1 

(1 + x 2 ) k 
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(x — 



Im(z) > 0. 


(35) 
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It easy to see that this integral converges for all k > 0. There are two cases 
to consider: the case where k > 0 is an integer, and the case when k > 0 is 
not an integer. The solution in the first case was derived independently by 
Chateau and Meyer- Vernet 2 and by Summers and Thorne. 1 The solution in 
the second case was derived by Mace and Hcllberg. 3 As will now be shown, the 
solution for arbitrary k > 0 is easily obtained by first finding the differential 
equation for I(z). This method of approach is an alternative to the method 
of contour integration used by the above authors. 

Differentiating equation (35) with respect to z yields 


n*) 



i 

(i + x 2 ) k 



(36) 


Multiplying this by z and then using the identity 


1 + 


2) 


X — Z 



(37) 


it follows that 


zl' + l 



1 

(i + x 2 ) k 



dx. 


(38) 


Differentiating this equation with respect to z, multiplying the result by z, 
and then using the identity (37) yields 


z 2 I" + 4 zl' + 21 



1 

(i + x 2 Y 


2x 2 

{x — zY 


dx. 


(39) 


Now differentiate (36) and add the result to the last equation to obtain 

r°° i 2 

(1 + z 2 )I" + 4 zl' + 21 = ■ dx. (40) 

J - oo (1 +X 2 Y \X-zY 

Integrating the remaining integral by parts and comparing the result with 
equation (38), the right hand side of (40) equals 

r°° I x 

- 2 < k - !) K dx = - 2 < k - W + ^ < 41 > 


Thus, the differential equation for I(z) is 

(1 + z 2 )I" + 2 (k + 1 )zF + 2 kI = 0. 


(42) 
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To obtain a unique solution, it is necessary to impose the initial conditions 


I{z = i) = /o,. 


I'(z = i) = 


2k tin 


k + 1 2 ’ 

where Jo is evaluated in Appendix A with the result 

: r( K + !) 


In ~ i V 5T 


Making the transformation 


t = 


t(k + i)' 

l + iz 


the differential equation (42) takes the form 

t{ 1 - t)J" + (k + 1)(1 - 2t)I’ - 2 kI = 0, 

where 


(43) 

(44) 

(45) 

(46) 

(47) 


I(t = 0) = (48) 

2 k 

J'(t = 0) = —J 0 . (49) 

k+1 

Equation (47) is now in the standard form of the well known hypergeometric 
equation. 

Recall that the solution of the hypergeometric equation 


z( 1 — z)F " + [c — (a + b + 1 )z]F' — abF = 0 (50) 

with initial conditions 

no) = i, (si) 

no 0) = y (52) 

is given by Gauss’ hypergeometric function F(a, b, c, z) = 2 Ei (a, b , c, z ). This 
is a single valued analytic function of 2 in the complex plane with a branch 
cut from 1 to 00 along the positive real axis. For fixed z it is analytic in a, 
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b, and c separately with the exception of simple poles at c = 0, —1, —2, ... In 
the unit disk it has the series representation 


ab z a(a + l)b(b+l)z 2 

F(a, b,c,z) = 1 + - - + ^ . 1 X 7 + 

c 1! c(c + 1) 2! 


\z\ < 1, 


(53) 


which also shows that it is symmetric in a and b. The series terminates if a 
or b G {0, —1, —2, . . . }, in which case F(a, b, c, z) is entire because it reduces 
to a polynomial in z. The properties of the function F have been thoroughly 
studied and documented. Therefore, for all intents and purposes it can be 
considered to be a closed form function in the same way that sin(z) and 
cos(z) are. 

Choosing 

a = 1, b — 2 k, c — k + 1, (54) 

equation (47) is equivalent to the hypergeometric equation (50); hence, the 
unique solution of (47) is given by 


I{t) = i v / 7T^7^yfyT 1 (l, 2k, k + 1, t). 

1 (k + 1) 

Using (46), a closed form expression for the integral (35) is 

- r (K+ s) „t, „ . . 1 + i.Z' 


(55) 


I(z) = I0 'ryrfy F ( 1 ’ 2k> K + 11 ~^r ) ■ 


(56) 


which is valid for all k > 0, in fact, for Re(«) > 0. Substituting this into 
(34), the plasma dispersion function is found to be 


Z«{£) 


»(«- I) 

k 3/2 


F( 1,2k,, k+ 1 ,t), 


(57) 


where 


l+iz 

2 



(58) 


It is important to note that the function Z K (£) inherits all the properties of 
F(a, b, c, z). Because F(a, b, c, z) has a branch cut from 1 to oo, the function 
Z k (£) is single valued analytic with a branch cut from ^ = —i\fk to £ = —ioo 
along the negative imaginary axis. Here, the £-plane is the complex frequency 
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plane, £ = is/ \k\v 0 , and the lower half plane consists of all the damped modes. 


The behavior in the neighborhood of the singularity £ = —i a/k can be 
obtained from the identity 7 


F(a, b, c, z) = (1 — z) c a b F(c — a,c — b,c, z). 


(59) 


This yields 


F[l,2K,K + l, 1+ n lZ ) = f 1 n 1 * ) f(k, 1 -K,K + 1, )■ (60) 


Using the formula 7 


it follows that 
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Therefore, one obtains the asymptotic relation 


(62) 
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6 Symmetry Relations and Taylor Expansion 

Using the well established theory of the hypergeometric function, the 
properties of Z K (£) may be derived with little effort. From the series (53), it 
is obvious that for a, b , and c real, 

F(a,b,c,z*)—[F(a,b,c,z)Y. (64) 

Therefore, replacing £ with — £* in equations (57) and (58), it follows that 

z.(-C) = -[z,K)]*. (65) 
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This relates the values of Z K (£) in the left half-plane to those in the right 
half-plane. The classical dispersion function for the Maxwell distribution has 
the same property. 

From the identity 8,9 


F[ 2a, 2b, a + b + l, 


i + z ^ r(a + b + |)r(|) 
^ r(a + i)T(6+|) 
r(a + 6+i)r(-|) 
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one has 
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+ 2kizF{1, k + — z 2 ). (67) 

Substituting the closed form expression' 

F(a, b, b , z) = F(b, a, b, z) — (1 — z)~ a , 
equation (67) takes the form 

1 + iZ\ r~_F(K + 1 ) 

7p 

2 ' 

Multiplying this equation by the leading coefficient in equation (57) yields 

iy/n r(/c) 
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zF{1,k + \,\,-z 2 ), (70) 


where z = f./y/F. Changing £ to — £ (or z to —z) and adding the result to 
the last equation yields the symmetry relation 


Z«{Z) + z K {-i) = 


2 iy/n r(/c) 

2 > 


K 1 / 2 F(k 


\D- 

K 


2 \ ~ K 


(71) 


This relates the values of Z K (£) in the upper half-plane to those in the lower 
half-plane. In the limit asK^oo one recovers the symmetry relation 

Z(f) + Z(-0 = 2 i-Cie-S 2 


(72) 
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for the Maxwellian dispersion function. 

Using the series expansion (53) for F(a, b, c, z ), equation (70) immediately 
yields the following series expansion valid for \z\ < 1: 




K 


2k- 1 

K 1 / 2 


zil + 


r(« — |) 

2k + 1 


l-z z ) + 


2k + 1\ (2k + 3 


5 


:-^ 2 ) 2 + --- 


(73) 


Here, the even powers are neatly summed into the first term. If necessary, the 
first term can be expanded using the binomial theorem. The series expansion 
can also be written in the more compact form 




. r(«) 

\ 1 / 2 r(K- \) 


(i + ^ 2 


— K 


A ^ T(k + w + 1/2) __ 2 

* 1/2r («-5)iS T(n + 3/2) 1 


however, the coefficients in (73) are much more explicit. 


z\<l, (74) 


7 Series Expansion for Large z 

Begin with the identity 7 

F(a, 6, c, z ) = — ^(-^)“ a F fa,a-c+l,a-b + l,~ 

1 (6)1 (c — a) \ z 

+ (V b - c + 1,6 - a + 1, i ] . 

I (a)T(c — 6) \ z 

Applying this to the last term in equation (69) yields 
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In addition, using the identity (68), this simplifies to 
1 + iz s 


F(l, 2k, k + 1, — 2^) = \/n 


T(k + 1) 2\— k 

= y(l + z 2 ) 


r(K + |) 
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Using the reflection formula 
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(77) 


(78) 


the last equation takes the form 
1 + iz 


f( i,2*,* + i,i±H) 
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+ y/ni tan(7 m ) ^ ^ \ ^ z(z 2 ) (k+1/2) ^1 + . (79) 

In order to get the phases right special care must be taken when combining 
the different power functions. This is because, for the principal branch of 
the power function, it is not always true that ( z a Y = z ai 3 or z±z% = (ziz 2 ) a . 
However, using the identities 

1 \ « 

(80) 


Z *+P = z a z /3 ? (]_ + z y* = z «(l + A) “ 

and 
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I —1 if Re(z) < 0 and arg(z) ^ +7T/2, 
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which are valid for all z 7^ 0, equation (79) can be written 
1 + iz y 
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2 
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It should be noted that for half-integer values of k the right-hand side is 
undefined. Otherwise, this equation holds for all z ^ 0 in the complex plane 
with a branch cut from —i to —ioo along the negative imaginary axis. For 
\z\ > 1, the hypergeo metric function on the right-hand side of (82) can be 
expanded using the series (53) to obtain the expansion 


F( 1,2«,«+1,1±^) = 

+ 'ita.n(TK)s g n[R e (z)]}(z 2 )~* (t + D 

+ +Ad (: l + 3-b- (4V) + 73 — vl + ' 

(k — 5 ) 2 1 j-«\ Z 2 / (|-k)(|-k)V Z 2 7 


(83) 


1 2 1 > 1. The factor sgn[Re(z)] ensures that the right-hand side is continuous 
throughout the upper half plane. This expansion can be arrived at inde- 
pendently using the theory of Laplace transforms as shown in Appendix B. 
Substituting (83) into (57) yields 


y/n T(k - \) 


1 + i tan(7T/c)sgn[Re(z)] >(z 


1 1 

~7k~z ) 1 + 


I -« 


+ 


z 2 ) (§-«)(§-«) 


1 + 4 

Z z 


+V + 


which is valid for \z\ > 1 and k ^ |, |, § , . . . It is important to note that this 
expansion is exact, not just asymptotic. 


8 Integration of the Differential Equation 

The differential equation (42) can be written in the symmetric form 

s{ (i+ +'~s[( i+ ++*>]H < 85 > 

This has the first integral 

h+*) l ~~£[h+*)’m\=c. 


( 86 ) 
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where 

/ °0 1 Y( k I l') 

(i + ^i ix = “Tw ' (87) 

Equations (86) and (87) can be verified by substituting (35) into (86), car- 
rying out the differentiation, and then using integration by parts and the 
identity (37). Therefore, the solution can be written 

I(z) = C K (l + z 2 y K f (1 + t 2 Y~ l dt, (88) 

J i 

where the lower limit of integration is chosen such that 

lim I(z) = * 2 | • (89) 

1 [K, + 1 ) 

If k = 1/2, the integral (88) has the closed form solution 


I(z) 

= -2(1 + z 2 y 1/2 sinh- 1 ^)-! , 

(90) 

where 

Sml1 1(2)= / (1 + (2)1/2 

(91) 

Equivalently, this can 

be written 

/ \ 



sinh _1 ( 2 )) = log + V z 2 + 1 j, 

(92) 


where log ( 2 ) is the principal branch of the logarithm. Substituting the result 
(88) into (34), the plasma dispersion function can be written 

z««) = -^(1 + *T“ [ (1 + tT' 1 it, (93) 

where z = This holds in both the upper and lower half planes and 

therefore represents the analytic continuation of Z K . 
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9 Solutions for Half-Integer Values: k — n+| 

Setting k — n + | in equation (69) yields 

f ( li2K>K + 1 , l ± H )^£|±| (1 + 22r -/ 2 

+ (2n+ l)izF(l,n + 1, §,-z 2 ). (94) 
By repeated application of the identity 7 
(b — a)(l — z)F(a,b,c,z) = (c—a)F(a — l,b,c,z) — (c—b)F(a,b—l,c,z) (95) 
one may prove the formula 

F(i,„ + 1, f, -f) = i V r(» + i-y) r (n +i) m 

V ’ ’ 2 ’ ; 2^ r(n+l) r«+ -mh 7 

m = 1 x 7 v ^ 7 

I 1 r(l) r(rz + ^) 3 2\ 

+ 2T(n + l) r(|) + ^ ( 96 ) 

where n = 1, 2, 3, . . . Inserting this formula into (94) and using the closed 
form expression' 


F(l, 1, —z 2 ) = -(l + z 2 )- 1/2 sin^(iz) 

xz 


one finds 


F( l,2*,«+l,i±“) 


r> + 1) 

r(n + 1) 


(l + ^ 2 )- n - 1/2 |"l + — sin _1 (iz) 

7 r 


Ib-tl) y- r <" + 1 ~ m | fa( i + ( 98) 

r (™ + !) “i r(n + | - m) 


where n — n + Thus, from (56) one obtains the closed form solutions 


I(z) = -2(1 + * 2 )-"-V 2 sinlT 1 ^) - y 


os) 

^ T(n +f-m) 



J.J. Podesta — Plasma Dispersion Function for the Kappa Distribution 18 


n — 1, 2, 3, . . . The first few solutions are as follows; for k = |: 


I(z) = —2(1 + z 2 )~ 3/2 


sinh 1 (z) — 


-2z(1 + z 2 )~ , 


( 100 ) 


for k — §: 


I{z) = -2(1 + z 2 )- 5/2 


sinh 1 (z) 


Til 


-z(l + z 2 )- 1 -2z(l + z 2 )- 2 , (101) 
o 


and for k = (y 


I(z) = —2(1 + z 2 )~ 7/2 


sinh 1 {z) — 


1 fi 4 

- -^(1 + z 2 )- 1 - -z{ 1 + z 2 )~ 2 - 2z(l + z 2 )- 3 . (102) 
15 3 

Multiplying equation (99) by the leading coefficient in equation (34), the 
plasma dispersion function is given by 


Z «K) = ! t + 


(n + I) 1 / 2 T(n) 1 ^ T(n + | - m ) ' 


-7=(l + ^ 

\/7T 


2\— n— 1/2 


sinh AM 


7T2 


, (103) 


where z = and n = 1,2,3,... A more explicit form is obtained by 

making the substitution given by equation (92). Similar expressions were 
first obtained by Summers and Thorne 10 using different methods. 


10 Solutions for Integer Values: hi = n 

If k — n is a positive integer, a closed form solution can be derived using 
the identity (61) to obtain a hypergeometric series which terminates; how- 
ever, this solution lacks mathematical symmetry. It is preferable to proceed 
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as follows. The hypergeometric series (53) can be written 

T(1 + m) T(2n + m) T(n + 1) t T ' 


F(l, 2n, n + 1, t) = ^ 


m = 0 


T(l) T(2n) T(n + l + m)m! 


T(n + l)y^ T(2n + m) 
r ( 2 «) “" 0 r (n+l+m) 

T(n + 1) ^ T (m + n) tm _ n 


E 


T(2n) ^ T (m + 1) 

m=n v ' 

Splitting the sum into two parts 


(104) 

(105) 

(106) 


F( 1, 2 n, n + 1, t) — 


T(n + 1) 
T(2n) 


y- r (m + n) tm _ n _y^T(m + n) 

f J m I * ^ 


y, m = 0 


m\ 


m = 0 


m\ 


the hrst part can be summed by the binomial theorem 
T(m + a) z r ‘ 


y. 

m — 0 


T(a) ml 


= (1 ~z)~ a , 0, -1, -2, 


(107) 

(108) 


to obtain 


F{ 1, 2n, n + 1, t) = 


T(n + 1) 
T(2n) 


n— 1 


r(«)r*(i - ()-” - £ r ( m + n) r- 


ra=0 


mi 


Substituting t = (1 + iz)/2 this yields 

1 + iz\ T(n + l)T(n) f f 1 + z 2 ^ ~ n 


(109) 


F^l,2n,n + 1, — ^) 


T(2n) 


T(m + n) ( 1 + iz 


n — 1 

V 

^ T(n)m! 

m=0 v ’ 


. ( 110 ) 


Using the duplication formula (B6) found in Appendix B to evaluate the 
coefficient on the right hand side, one obtains 


F(l, 2 n, n + 1, 


1 + fa ' = 2^ T JKDL{ (1 + zV 


r(n + 1) 

n— 1 

E 


1 r(m + n) ( 1 + iz 
4 n 2-*/ r (n)m\ 

m = 0 v 


■ (HI) 
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Comparing this equation term by term with the symmetry relation [see equa- 
tion (69)] 


f(i, 2k, k + 1, + F (l, « + i, 2 ^ ) 

it follows that 


2^ 


r(/c + i) 

r(« + |) 


(i + ^ 2 )- K 

(H2) 


F^l, 2 n, n + 1, 


1 = i / 2 rv + i ) 1 

2 ) r [n + \) 4 n 


^-4 T(m + n) 
T (n)m\ 

m = 0 v ’ 


f l + iz 

V ~2“ 


or, equivalently, 


m—n 


(113) 


f(l,2n,n + 1, ' if ') 


7T 


1/2 


n 


2 2n 1 Tin + 1) ' ml 

v m=0 


n— 1 

vE 


r (m + n) / 1 — iz 


Substituting this into equation (56) yields 

n— 1 


'M = £rE 


m= 0 


n + m — 1 \ ( 1 — iz 


m 


(114) 


(115) 


n = 1, 2, 3, . . . This same result is obtained by using the theory of residues 
to evaluate the integral (35) along a large semicircle in the lower half plane. 2 
Because (115) is a rational function of z, it immediately gives the analytic 
continuation of I(z) throughout the complex plane. 

Using equation (110) to rewrite (113) in the form 


’(l,2n, n+1, 


1 + iz\ T(n + l)T(n) ^-4 T(m + n) (1 — iz 


T(2n) 


E 


— ' Tin)ml 

m = 0 v ' 


(116) 


it follows from equation (57) that the plasma dispersion function is given by 


7 - - ( n ~ I) r W r (' m + n ) ( 1 ~iz \ m U 

* n 1 / 2 T(2n) ^ o ml \ 2 J 


(117) 


where z = £/a Jn and n — 1, 2, 3 . . . 
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11 Equivalent Integral Representations 


Consider the integral defined for Im(^) > 0 by 

r°° i l 

m = 


L 00 (1 + x 2 ) k (x-z) 


dx 


(118) 


and by analytic continuation for Im(z) < 0. By means of some elementary 
changes of variables, one may derive the following equivalent representations: 


I{z) = 

I(z) = 


oo(l + a;2 ) K (x 2 — z 2 ) 

r i i 


dx, 


dx, 


l-oo (cosh a;) 2 * 1 [sinh(x) z] 

l ^ 

J (-)=/ • r ■ , 2^ 


I(z) = 


Loo (c° s hx) 2K 1 [sinh 2 (x) — z 2 \ 

dx, 


' — OO 

/»oo 


2 \ « 

£ \ Z 


1 + X 2 


1 — Z 2 X 2 


, i . , , 9r z cosh x 

I{z) = / (tanhx) ^ — dx, 

J _ oo 1 — z 2 sinh x 

/»oo 

I(z) = z / t K_1 / 2 (l -f t) _K (l — z 2 t ) -1 dt, 

Jo 

Thus, by equation (55), all of these integrals are equal to 

T , . . / —T(k + 1J2) ( 1 + iz 

I(z) = — -F ( 1 , 2k, k + 1 , 

V ; V T(/c + l) V 2 


(119) 

( 120 ) 
( 121 ) 
( 122 ) 

(123) 

(124) 

(125) 


In addition, I(z) can be represented as a Laplace transform. If z — is 
and Re(s) > 0, then one has the identity 


= i [ e~^ s+ix)t dt. 

X — Z S + IX J o 


(126) 


Substituting this into (118) and changing the order of integration yields 


where 


/(f) = 2 


f{t)e st dt, 

(127) 

cos(xt) 

7 - 77 - dx. 

(1 + x 2 ) k 

(128) 
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Using the well known integral representation of the modified Bessel function 
of the second kind: 


K u (t) = 


r(z/ + 1/2)2" 


cos (xt) 


'k x !H v J 0 (1 + x 2 Y +1 / 2 


dx, t > 0, (129) 


it follows that 


= i>0 ' 


This establishes the identity 

I(z — is) = i2zf[/(t)], Re(s) > 0, 


(130) 


(131) 


where /(f) is defined by the previous equation and SP is the Laplace transform 
operator. 


12 Conclusions 

Starting with the definition of the plasma dispersion function in terms 
of the integral (35), a second order differential equation is derived that im- 
mediately allows the dispersion function to be expressed in terms of Gauss’ 
hypergeometric function. Using the well known properties of the hypergeo- 
metric function, the many properties of the plasma dispersion function are 
then derived in a unified manner. A summary of the results is provided in 
Section 2. 


Appendix A: Evaluation of the Integral Jo 

By definition 


J 0 = I{z = i)= f 


dx. 


l-oo (1 + X 2 Y (x-i) 
Multiplying the integrand by (x + i)/(x + i), this becomes 

1 


Jn — 2 i 


'o (1 + X 2 Y +l 


dx, 


(Al) 


(A2) 
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and with the change of variable t = x 2 one obtains 


In = i 


t-y 2 


./u (l+t) K+1 

The remaining integral is the Beta function 


dx. 


B(p, q) = 


ro ° t?- 1 r(p)r( q ) 


I o (1 + t)P+<i^ r {p + qY 

with p — 1/2 and q — k + 1/2. Therefore, it follows that 

r(/c + 1/2) 


In = W7T 


T(/c + l) 


The calculation of /'(z = i) is similar. 


(A3) 


(A4) 


(AS) 


Appendix B: Alternative Derivation of the Ex- 
pansion for \z\ > 1 


The purpose of this appendix is to give an independent derivation of the 
expansion (83) for \z\ > 1. The starting point is the representation of that 
function as a Laplace transform. Using the definition 


K v {x) 


TT I- V {x) ~ I„(x) 
2 sin(7rz/) 


(Bl) 


together with the series expansion for the modified Bessel function I v (x), one 
obtains 



71 


E 


2sin(7rz/) I -4-^ T(n — v + l)n! \2 


x 


2 n 


E 


^ Tin + v + l)n! 
n= o v ' 


2n+2tA 

d } 


(B2) 



J.J. Podesta — Plasma Dispersion Function for the Kappa Distribution 24 


where x > 0 and v > 0 is not an integer. By a well known theorem, 11 the 
Laplace transform of this function has the asymptotic expansion 


& 



7T 



2 sin(7rz/) 



r(2n + l) 1 1 

r(n — v + 1) n!2 2n s 2n+1 


OO 


-£ 

n = 0 


r(2n + 2z/ + l) 1 1 

r(n + v + 1) n\ 2 2n + 2u s 2n+2v+l 


(B3) 


where s — > oo along any ray e ld with \9\ < ir/2. In fact, the following analysis 
shows that this series converges in the right half-plane whenever |s| > 1 and, 
therefore, equality holds throughout this region. However, this fact is not 
needed now. Setting v = k— 1/2, then from the relations (131) and (56) one 
obtains 


F[ 1 , 2k, k + 1 , 


2 

71 


1 F( ' t+ h/(«o 


i^/nY(K + \)' 


K 


IT E 


r(2n + 1) 1 


2sin[(/c — \)ir\ T(k + \) [ “ T(n — k + |) n!2 2n s 2n+1 


-E 


T(2n + 2k) 


Yin + k + i) n!2 2 "+ 2 «- 1 s 2n + 2ft 

n=0 v 2 ' 


. (B4) 


where s — » oo along any ray e* e with |d| < 7 t/ 2. It will now be shown that 
this agrees with the expansion (83), that is: 


F[ 1 , 2k, k + 1 , 


1 — s 


2 

t 1 + i tan ( 7rK ) s § n [ Re (* s )]] (- s2 ) _K ( 1 


+ 


K 


[K - 2 ) S 


177 1 1 + 


1 

WT I + 


1 3 

2 ‘ 2 


« (§-»)(!«») V* 


The even and odd power terms are compared separately, 
terms. Using the duplication formula 



First the even 


r(2 z) = !r 1/2 2 2 '- 1 ro)r (z + 1 


(B6) 
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and the binomial theorem, the even terms in (B4) can be summed as follows 


7T k 1 T(2n + 2n) 1 

sin[(/c — |)7r] r(/c + |) n\ T(n + k + |) 2 2n+2K_1 s 2n+2K 

y/F /cr(/c) ^^T(n + K) 1 

cos(k7t) V(k + ^ r(/c)n! s 2n+2K 

v ' n = 0 

= « r (*0 c - 2 k ( 1 _ 

cos(k7t) T(k + \) \ s 2 ) 


— AC 

. (B7) 


This will now be compared with the even terms in (B5): 


\/^ r(^+ 1) t 1 + i tan ( 7rK ) s § n [ Re (* s )]] (-s 2 ) 

It is easy to show that 

(~s 2 )~ K _ ( e +i7rK if 0 < arg(s) < vr/2 

s~ 2k \ e _l7rK if —n/2 < arg(s) < 0, 

therefore, 

s~ 2k 

7 — - = |"l + i tan(7r«;)sgn[Re(is)]l (— s 2 ) _K , 

cos(k7t) l j 

which establishes the equality of (B7) and (B8), that is, for the even power 
terms in (B4) and (B5). Similarly, using the duplication formula (B6) the 
odd terms in (B4) are easily shown to be equal to the odd terms in (B5). 
This gives an independent proof of (83) in the upper half plane which, by 
means of the symmetry relations, extends to the whole plane. 


(B8) 

(B9) 

(BIO) 
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